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Color Atomic Force Microscopy with on-the-fly Morse parameters mapping 
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Atomic Force Microscopy has enabled 2D imaging at the sub-molecular level, and 3D mapping 
of the potential field. However, fast identification of the surface still remains a challenging topic. 
In this paper, as a step towards implementation of such function, we introduce a control scheme 
and mathematical treatment of the acquired data that enable retrieval of essential information 
characterizing the potential field, leading to fast acquisition of images with chemical contrast. The 
control scheme is based on tip sample distance modulation at an angular frequency of uj, and null 
control of the lu component of the self excitation frequency of the oscillator. It is demonstrated 
that the control is robust in UHV for a frequency well as small as a few Hz/MHz, and that the 
mathematical treatment results in satisfactory identification of the potential field. Morse potential is 
chosen as a case study for identifying the Morse parameters per pixel. Atomic features with similar 
topography were distinguished by differences in the parameters. The decay length parameter was 
resolved with a resolution of 10 pm. The method was demonstrated on quenched Silicon at a scan 
rate comparable to normal imaging. 

PACS numbers: 37,68,73,82 


In dynamic mode Atomic Force Microscopy, an oscilla¬ 
tor, such as a cantilever or a tuning fork is used to vibrate 
the probing tip in proximity of the sample surface. There 
are basically two types of operational modes of imaging, 
which are: (i) the constant frequency shift mode [1], and 
(ii) the constant height mode [2, 3]. The former, although 
widely used for imaging, risks loosing a large portion of 
information, since the obtained equi-frequency shift sur¬ 
face is dependent on the imaging parameters and do not 
necessarily retain information of the potential landscape. 
The latter has the possibility of yielding more informa¬ 
tion, since the logged frequency shift can be used to cal¬ 
culate the force or potential [4, 5]. The mode, ideal as it 
may sound, can only be used in low drift conditions, due 
to the absence of tip sample distance regulation, and is 
often time consuming. In the case of CO terminated tips, 
it also has the difficulty in probing the surroundings of 
atomic or molecular protrusions due to masking effects 
and possible tip tilting mechanisms [6, 7]. 

With a view to solving the problems listed above, and 
to implement microscopy with chemical contrast, we in¬ 
troduce here a control scheme that sets the working point 
to the vicinity of the bottom of the frequency curve (FC), 
and mathematical treatment to parametrize the acquired 
data as parameters that define the potential field between 
the tip and the sample. The control of tip sample dis¬ 
tance is accomplished by applying position modulation 
of the sample in the tip sample direction with an angular 
frequency of cj. Tip sample distance is regulated to null 
the luj component of the oscillation. As a result, in the 
case of low dither amplitudes, the bottom of the FC is 
maintained as the working point. For larger dither am¬ 


plitudes, the working point moves away from the sample 
surface due to asymmetry of the FC. 

FC is uniquely defined per site for a given amplitude 
of drive A. Calculation based on the acquired data is 
facilitated by the boundary conditions assured through 
the null control of Afi^. Such position modulation tech¬ 
nique was first introduced to SPM as a differential imag¬ 
ing technique [8, 9] then applied to lateral position con¬ 
trol [10-12] or tip sample distance regulation [13, 14]. 
Tip sample distance modulation was also used for syn¬ 
chronous detection of the harmonics of modulated fre¬ 
quency shifts in order to reconstruct full frequency shift 
curves between the tip and the sample [15]. Here, po¬ 
sition modulation of the tip sample distance is set to a 
few angstroms with a view to scanning, in the vertical 
direction, the skin depth of the short range forces. 

Ideally, not to loose any information, all the higher 
harmonics of the oscillation should be logged and used 
for the characterization of the surface. In actual imag¬ 
ing, depending on the amplitude of drive of the funda¬ 
mental mode of oscillation, the higher harmonics may 
have much smaller amplitudes, rendering measurement 
with reasonable frequency noise difficult. Here as a case 
study, we adopt the Morse potential, and identify the 
Morse parameters for every pixel of imaging by logging 
the DC and the 2uj component of the modulated oscilla¬ 
tion. Suitability of working with a limited order of higher 
harmonics needs to be evaluated through the fitting abil¬ 
ity of the identified parameters to actual measurement, 
and through imaging of various samples. 

From its easy physical grasp, with its three fitting pa¬ 
rameters, the Morse potential is a widely used model to 
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describe inter-atomic potentials [16, 17], gases and sur¬ 
faces [18-21], and also widely used for inter-atomic in¬ 
teractions in molecular calculations [22-28]. The Morse 
potential is given by: 


2e--(1) 

Where Ei^ is the energy depth of the bond, Zq is the 
equilibrium distance, and A is the decay length that 
parametrizes the width of the well. These three param¬ 
eters serve as chemical fingerprints of the atoms on the 
surface. 

The frequency shift of a cantilever that oscillates with 
amplitude A is given by [29-31]: 
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/o and k are eigenfrequency and the stiffness of the 
flexural mode. Modified Bessel functions of the first kind 
In{x) are used in preference to Kummer functions for 
ease of manipulation of the equations. For the specific 
case where ^4 << A, this frequency shift is proportional 
to the second derivative of the potential [29]. In the 
general case of amplitude of ^4, if a dither amplitude D 
is applied around the position Zj:fQ at angular frequency 
uj to keep the condition A/i^ = 0, we obtain: 
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Solution of Eq. (3) is solved for the measured center 
of dither Zdo 


With 
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Figure 1: Schematic of the parameter mapping method. 

LDV: Laser Doppler Vibrometer. PI: Proportional Integral 
Feedback controller. 
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/, gi and g 2 are dimensionless functions of a and p. 
In the range where A/dc and Af 2 uj are monotonic, (i.e. 
for every a and for P <« 2), the system can be uniquely 
solved for and A. A is injected into Eq. (4) to obtain 
Zq. There are no general analytical expressions for the 
solution of the system but it can be solved numerically 
[32]. Eor ^4 << 1 we have : 


hl{a) and h2{p) (with a = A/X and p = D/X) are 
functions of A and D and therefore dimensionless func¬ 
tions indicating the relative position of Zjjq with respect 
to Zq. Zdo is the controlled value of the tip-sample dis¬ 
tance, and can be seen as a universal topography for it 
is expressed uniquely as a function of the Morse param¬ 
eters. 

Since Af is a periodic signal around the working point, 
the DC and 2uj components of Af are used to derive the 
Morse parameters: 


Nmin - J ( 8 ) 

and, 

F = — (9) 

giving a rough idea on the physical meaning of local 
minima of frequency and force. 

Eigure 1 depicts the schematic of the AEM. The A/i^ 
signal is fed to the z feedback control loop to keep its 
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Figure 2: Determination of Eh and A versus dither ampli¬ 
tude on top of Si(lll) from experimentally obtained frequency 
components A foe, Af 2 uj and solving equations (6b) and (6d). 
(a) Experimental values of A foe and Af 2 uj- (b) Effective 
bonding energy Eh. (c) Effective decay length A. Acquisition 
parameters were f 2 nd ~ 1-99 MHz, fend = 2,400 A 

= 7 A, uj = 1 kHz. 


value close to zero. Angular frequency uj is set below the 
frequency of self-excitation around 2 MHz and above the 
cut-off frequency of the feedback control. 

Figure 2 shows dependence of calculated E^ and A on 
D while executing null control of A/i^. The values are 
more or less insensitive to which is a good sign, since 
the identified Morse parameters should be intrinsic to the 
site and not affected by experimental conditions such as 
D. Small D should be avoided since the tip apex does not 
make sufficient excursion within the short range force, 
making the 2uj signal too weak. Large D should also be 
avoided since Afoc converges to zero due to excessive 
averaging of the FC well. 

Figure 2(a) depicts the measured values of Afjjc and 
Af 2 uj- In simulation and in actual measurement, signals 
/A/dci A/ 2 u; cross around 0.75A, which gives a guideline 
on the favorable value of D for minimizing the measure¬ 
ment noise. The values in Figs. 2(b) and 2(c) are snap 
shots of the calculated Morse parameters. The values 
appear to be scattered for amplitudes of D far from the 
cross-over point but otherwise give compatible noise level 
for scanning experiments. Imaging shows that the values 
are site dependent and sufficiently low-noise to differen¬ 
tiate the atomic features. 

Figures 3 and 4 depict signals A/dc^ Af 2 uj and Zdo 
acquired on quenched Si(lll) while keeping the luj com¬ 
ponent of Af to zero. A and Zq were numerically 
calculated from the acquired data. Details on the AFM 
can be found in the literature [33-35]. 


Figure 5 depicts a measured FC and a FC generated 
from a set of Morse Parameters chosen to fit the measure¬ 
ment. The fit is satisfactory around the bottom of the 
FC where the working point is modulated by but devi¬ 
ates further away from the sample surface. This may be 
due to the commonly known fact that the Morse Poten¬ 
tial is better adapted for smaller inter-atomic distances 
[17]. Adopting a potential field model with more param¬ 
eters may result in a better curve fit of the FC. What we 
propose is to apply the method to various samples to be¬ 
gin with, and if it proves to be limited, move onto other 
potential models using measured higher harmonics. One 
may note that the fit is already satisfactory within the 
range of the position modulation and that larger D 
may be needed to identify parameters that fit the entire 
FC. 

From the experimental results, we have demonstrated 
the following: (i) atomic resolution is obtained in UHV 
for the proposed control method, and control is stable 
even for a small local minima of frequency of 5 Hz over 2 
MHz. (ii) Calculation from the measured data to the 
Morse parameters took 1 sec per pixel, but could be 
shortened to 2 ms by preparing a look-up table for the 
conversion, (iii) values of E^ and A tend to be larger than 
expected for a Si-Si pairwise potential, (iv) if atomic fea¬ 
tures take the same colour, it means that the set of Morse 
parameters are the same. However, care is needed in in¬ 
terpretation of the images. Atomic features A and B are 
similar atomic features in the topography Zdo image, 
but different in the final colour image due to difference 
of A. Difference in A could be resolved in the 10 pm or¬ 
der. Features B and C have more or less the same A and 
Eh but differ in colour due to the topographic difference. 
The derived Morse parameter Zq is directly affected by 
topography such as atomic steps. The sample should be 
atomically flat, or the topography channel (Red) should 
be used with care. The method proposed here starts off 
with measured FC and derives a set of Morse parame¬ 
ters that give a satisfactory curve fit between measured 
and calculated FC. The measured FC is a result of var¬ 
ious interaction between the tip and the sample, so due 
to such effects as back ground forces and relaxation, the 
derived Morse parameters may deviate from the pairwise 
potential expected of the tip apex atom and the sample 
atom. Deviation needs to be assessed in a quantitative 
manner to deduce Morse parameters close to the pairwise 
potential. If we limit the discussion to an atomically fiat 
surface, the discussion can be simplified to the following: 
(i) tip volume acts to increase van der Waals force, result¬ 
ing in larger Eh and larger relaxation of the tip and the 
sample. The latter results in a larger A due to fishing-out 
and pushing-in of the sample atom. Through calculation, 
placing a molecule at the apex of an AFM tip has the fol¬ 
lowing effects in view of the method proposed: (i) Due to 
van der Waals attraction, effective Eh is minimal when 
tip radius R is zero, and increases with R. Effective Eh is 
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Figure 3: Morse parameters mapping on quenched Si(lll). (a) Measured quantities using Eqs.(4), (6b) and (6d). (b) 

Calculated Morse parameters, (c) RGB Morse parameters mapping. Image size is 2.2 x 2.2 nm^. Acquisition parameters were 
f 2 nd ~ 1-99 MHz, fend = 2,400 A = 7 A, u = 1 kHz and D = 4 A. 
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Figure 4: Experimental parameter mapping results on quenched Si(lll). The dimension of each picture is 6x5 nm^. Scale units 
are stated above the scale bars. Measured values : (a) A foe, (b) Afeo; and (c) Zdo- Calculated Effective Morse Parameters : 
(d) Efe, (e) A, (f) Zq and (g) RGB combined picture (R = Zo; G = A; B = Eb) with distinguishable atoms, (h) reconstructed 
Morse potential profile for atoms A, B and C. Acquisition parameters were / 2 nd = 1.991 961 MHz, fend = 2,400 A = 

7 A, uj = 1 kHz and D= 4 A. 


2 times larger for R^IO nm [32]. Functionalising the tip 
apex with a molecule [36] greatly lessens the increase of 
effective due to R, (ii) for a rigid body model A is much 
less sensitive to increase of R. However, if relaxation of 
the sample and the tip are taken into consideration, a 
sharper tip or a molecular functionalized tip is effective 
in lessening the difference between effective A and that of 
the pairwise potential. As a rough calculation, assuming 
the spring constant of the sample atom to its surrounding 
to be 10 N/m, attractive force of 160 pN results in 16 pm 
of relaxation that is not accounted for in calculating A. 
Calibration of the same order is necessary in the repul¬ 
sive regime. It was demonstrated that the method allows 


resolving the difference in effective A in the order of 10 
pm. This difference can be ascribed to at least two ef¬ 
fects, namely, (i) the difference of the pairwise potential, 
and (ii) the difference of bonding rigidity of the sample 
atom to its surroundings. Further study and simulation 
is needed to access the ratio of the two contributions to 
the contrast mechanism. As mentioned above, introduc¬ 
ing a molecular terminated tip, such as an adamantine 
molecule terminated tip [37], will act in favour in lessen¬ 
ing the difference between effective E^, A and the pairwise 
Eb and A. Using a well defined molecular tip will improve 
the repeatability of the measurements, and quantitative 
evaluation of the sample. 
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Figure 5: Curve fitting between measured FC and FC gen¬ 
erated from effective Morse Parameters, Eb = 2.2 eV and 
A = 4.6 A and using Eq. (2). Same paremeters were used as 
Fig. 2 

The proposed method is robust, and in many cases, 
tip change does not occur during acquisition of an im¬ 
age. This makes relative comparison of the atomic fea¬ 
tures within the image possible. However, to establish 
quantitative observation of the surface, a reference tip 
termination method needs to be sought for, and the de¬ 
rived physical values studied in terms of resolution and 
absolute values. Effects of combination of neighbouring 
atoms and subsurface compositions need also be studied, 
and a database constructed for a better understanding. 
The relatively simple and robust method introduced here 
offers such possibilities. We hope that the method will 
serve as a milestone towards such efforts. 
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